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Introduction 


This  report  documents  a  set  of  FORTRAN  programs  and  subroutines  for 
analyzing  stratigraphic  time  series  data.  These  programs  were 
written  as  part  of  the  author's  Ph.D.  dissertation  (Schiffelbeln, 
1984)  and  can  be  considered  as  an  appendix  to  that  work.  This 
documentation  is  divided  into  two  sections.  The  first  section 
follows  the  organization  of  the  author's  thesis.  Each  project  is 
briefly  explained  to  provide  a  context  for  the  programs.  All 
algorithms  used  in  the  project  are  also  explained.  The  second 
section  provides  an  alphabetical  listing  of  the  algorithms  referred 
to  in  the  first  section.  All  programs  were  written  for  the  UCSD 
Vax  11/780  (VMS)  computer.  Subroutines  referred  to  in  the  programs 
but  not  listed  (CARTESIAN,  CLEAR,  COMPLEX,  CONVOLVE,  FFT, 
FFTINVERSE,  INTERPOLATE,  POLAR,  POWER,  POWERDB,  PROLATE,  REALPART, 
ZEROMEAN)  are  from  a  library  kindly  made  available  to  the  author  by 
Alan  Chave  and  are  described  and  listed  in  Chave  (1980). 
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The  Effect  of  Benthic  Mixing  on  the  Information  Content  of 
Deep-Sea  Stratigraphic  Signals 

Abstract  (Chapter  3  of  the  dissertation): 

Benthic  mixing  or  bioturbation  affects  sediments  in  a 
number  of  ways,  including  1)  the  production  of  trace  fossils, 

2)  mechanical  and/or  chemical  alteration  of  the  sediment,  and 

3)  filtering  or  smearing  stratigraphic  signals.  Since  mixing 
alters  both  the  slopes  and  amplitudes  of  any  recorded  events, 
some  knowledge  of  the  process  is  essential  for  correct 
interpretation  of  the  signals.  Particularly  in  the  light  of 
recent  trends  towards  high-resolution  stratigraphy  and  signal 
unmixing,  the  frequency  characteristics  of  the  benthic  mixing 
filter  must  be  understood  to  know  which  types  of  signals  can 
be  detected  after  mixing.  Ash  and  tektite  profiles  in 
deep-sea  cores  from  various  geographic  regions  are  examined  in 
the  frequency  domain  to  determine  the  characteristics  of  the 
benthic  mixing  process.  Analysis  of  these  profiles  Indicates 
that  even  signals  from  cores  having  sedimentation  rates  as 
high  as  7  cm/ka  will  show  severe  attenuation  of  frequencies 
higher  than  0.35  cycles/ka  (periods  shorter  than  2.9  ka) 
resulting  in  loss  of  the  ability  to  resolve  closely-spaced 
events.  Most  signals  will  experience  much  more  serious 
high-frequency  attenuation,  however.  The  severity  of  high 
frequency  loss  is  directly  related  to  sedimentation  rate  (R  « 
0.97  for  the  cores  examined),  suggesting  that  this  is  the  most 


important  variable  when  considering  a  core  for 
paleoceanographic  or  paleoclimatic  study. 
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Program  PROJECT:  This  program  computes  the  power  spectrum  of  an 
impulse  response  function  (IRF).  The  IRF  is  stored  in  some  input 
file.  The  power  spectrum  is  given  in  decibels  and  stored  in  a 
user-specified  output  file. 

Program  GSTEST  is  a  simple  driver  program  for  the  subroutine  GSIRF. 
The  user  is  asked  for  the  benthic  mixing  parameters,  which  are 
passed  to  GSIRF,  and  the  mixing  function  is  generated  and  stored  in 
a  user-specified  file. 

Subroutine  GSIRF  uses  the  approximation  equation  of  Officer  and 
Lynch  (1983)  as  a  solution  to  the  benthic  mixing  equation  of 
Guinasso  and  Schink  (1975).  The  subroutine  is  accessed  by  a 
statement: 

CALL  GSIRF  [generated  IRF  vector,  mixed  layer  thickness 

2 

(cm),  diffusion  coefficient  (cm  / ka) ,  sedimentation  rate 
(cm/ka)].  The  output  is  the  benthic  mixing  impulse 
response  function. 


B.  Extracting  the  Benthic  Mixing  Impulse  Response  Function:  a 
Constrained  Deconvolution  Technique 

Abstract  (Chapter  4  of  the  dissertation) : 

Benthic  mixing  or  bioturbation  acts  as  a  low-pass  filter 
on  stratigraphic  signals,  altering  both  the  apparent  rates  and 
amplitudes  of  recorded  events.  The  mixing  process  is 
represented  as  a  time-invariant  convolution,  and  uses  an 
impulse  response  function  (IRF)  parameterized  in  terms  of 
mixed  layer  thickness  and  mixing  intensity  (diffusion 
coefficient).  Such  parameterization  facilitates  construction 
of  the  inverse  (deconvolution)  model,  which  is  used  to  examine 
the  effects  of  incorrect  IRF  characterization  on  deconvolution 
results.  Results  show  that  mixing  intensity  is  as  important  a 
parameter  as  mixed  layer  thickness  when  unmixing  a  strati¬ 
graphic  signal. 

Stable  oxygen  isotope  records  from  two  different  species 
of  planktonic  foraminifers  with  changing  relative  abundances 
in  the  same  core  will  show  apparent  offsets  in  the  timing  of 
events  due  to  mixing  (Hutson,  1980).  The  offset  can  be  used, 
with  the  appropriate  unmixing  equation,  to  constrain  mixed 
layer  thickness  and  mixing  intensity  based  on  analysis  of 
actual  isotope  signals.  The  assumption  is  that  any  offsets  in 
the  records  of  the  two  species  are  totally  a  product  of  the 
mixing  process. 

The  technique  was  applied  to  a  high-resolution  glacial- 
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interglacial  record  (Termination  II)  from  an  equatorial 

Pacific  piston  core.  Multiple  lateral  subsamples  were  stacked 

to  increase  signal  to  noise  ratio.  A  multidimensional 

nonlinear  optimization  routine  (SIMPLEX)  was  used  to  minimize 

an  error  function  related  to  stratigraphic  offsets  between  the 

two  species.  A  solution  to  this  problem  yielded  a  mixed  layer 

2 

5  cm  thick  and  a  mixing  intensity  of  7  cm  / ka.  This  mixed 

14 

layer  thickness  is  at  the  low  end  of  estimates  based  on  C 

profiles  in  box  cores  from  the  same  geographic  area  (Berger 

and  Killingley,  1982).  The  mixing  operator  is  probably 

conservative,  and  does  not  take  into  account  that 

sedimentation  rates  were  higher  during  the  last  glacial  period 

than  during  the  Holocene.  This  may  explain  the  inability  of 

the  algorithm  to  completely  remove  the  stratigraphic  offset 

between  the  two  isotope  records.  The  unmixed  record  (using 

2 

the  parameters  5  cm  and  7  cm  /ka)  shows  a  deglaciation 
overshoot  which  is  distinctly  smaller  than  the  "meltwater 
spike"  proposed  by  Berger  et  al.  (1977),  i.e.,  less  than 
0.1  °/oo  versus  a  range  of  0.3  °/oo  to  1  °/oo. 


Program  HUTSON:  This  program  performs  a  benthic  mixing 
(convolution)  operation  on  an  input  stable  isotope  signal. 

The  mixing  takes  into  consideration  the  change  in  relative 
abundance  of  the  isotope  carrier  (foraminifer  species)  as 
described  in  Hutson  (1980).  The  user  is  asked  for  input  file¬ 
names  for  isotopes,  abundances  and  the  mixing  operator.  The 
mixed  record  is  stored  in  a  user-specified  file. 

Program  DECOTEST:  This  program  is  a  simple  driver  for 
time-domain  deconvolution  of  a  stable  isotope  signal.  The 
user  is  asked  for  the  isotope  record  infile,  the  mixing 
function  infile,  and  a  constant  that  controls  the  sensitivity 
of  the  unmixing.  Values  between  0.001  and  0.05  were  found 
appropriate  for  most  real  and  simulated  isotope  signals,  with 
larger  numbers  giving  less  sensitivity.  The  program  accesses 
subroutine  DECO,  which  performs  the  deconvolution.  The 
unmixed  record  is  stored  in  a  user-specified  file. 

Program  HUTDEC:  This  program  is  the  inverse  of  program 
HUTSON.  Like  DECOTEST,  this  routine  deconvolves  an  isotope 
signal,  but  the  relative  abundance  of  the  signal  carrier  is 
also  considered  in  HUTDEC.  The  program  operates  as  described 
for  DECOTEST  above. 

Program  OFFSET:  This  is  a  large  and  relatively  complicated 
driver  program  that  attempts  to  remove  stratigraphic  offset 
between  two  stable  isotope  signals  from  different  foraminifer 
species  in  a  single  core  by  deconvolving  the  signals  with 
various  unmixing  functions.  (Refer  to  the  thesis  for 
details.)  The  heart  of  the  driver  is  a  simplex  optimization 
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algorithm  (Nelder  and  Mead,  1965)  modified  from  Daniels 
(1978).  The  fit  of  the  two  unmixed  curves  is  optimized  by 
adjusting  the  benthic  mixing  parameters.  The  program  examines 
stratigraphic  offset  by  calling  subroutine  ERROR,  which,  by 
calling  GSIRF  and  DECO,  unmixes  the  two  signals  and  examines 
closeness  of  fit. 

Subroutine  ERROR:  This  algorithm  calculates  a  value  for  the 
optimization  error  function  as  described  for  OFFSET.  This  is 
a  specialized  subroutine  which  accepts  two  isotope  signals 
(different  foraminifer  species  from  the  same  core)  their 
relative  abundances,  and  the  benthic  mixing  parameters.  Four 
calls  to  DECO  unmixes  the  records  and  calculates  the  residual 
between  the  two  unmixed  signals  (error  function). 

Subroutine  DECO:  This  algorithm  performs  a  time-domain 
deconvolution  of  a  measured  or  theoretical  stable  isotope 
signal.  The  subroutine  is  accessed  by  a  statement . 

CALI  DECO  [IRF  vector,  number  of  points  in 
input,  input  vector,  output  vector,  decon¬ 
volution  sensitivity]. 

DECOTEST  is  a  simple  driver  program  for  using  this  subroutine. 
Subroutine  CNVLV:  This  subroutine  convolves  two  vectors  and 
stores  the  result  in  a  third.  The  input  arrays  have  lengths  N 
and  M,  with  the  output  of  length  N  ♦  M  -  1. 

Subroutine  ZERO:  This  algorithm  sets  the  mean  value  of  a 
vector  to  zero.  The  original  mean  value  is  saved  as  AMEAN. 
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Bruhnes  Dissolution  Cycle:  Effects  on  Oxygen  Isotopes 
Sedimentation  Rates,  and  Signal  Spectra 


Abstract  (Chapter  5  of  the  dissertation): 

A  long  wavelength  cycle  (400-500  ka)  has  been  found  in 
stable  oxygen  isotope  records  from  Pleistocene  deep-sea 
sediments  in  the  Pacific  Ocean.  This  cycle  shows  an  Increase 
in  amplitude  with  increasing  water  depth,  correlates  in  both 
phase  and  duration  with  the  Brunhes  dissolution  cycle  of 
Adelseck  (1980)  and  is  apparently  a  result  of  differential 
dissolution.  This  dissolution  cycle  has  significantly 
affected  sedimentation  rates,  particularly  in  deep  water 
sediments.  Sediment  loss  through  dissolution  can  be 
quantified  by  comparing  isotope  stage  lengths  with  some 
shallow  reference  core  or  by  using  %CaC0^  data.  Examination 
of  a  number  of  equatorial  Pacific  cores  suggests  that  even 
sediments  recovered  from  a  water  depth  of  3500  m  have 
undergone  significant  distortion  as  a  result  of  this  disso¬ 
lution.  Long-wavelength  sedimentation  rate  perturbation 
results  in  a  widening  of  spectral  peaks,  with  a  consequent 
decrease  of  spectral  resolution,  if  some  depth  scale 
correction  is  not  made. 

The  Brunhes  dissolution  cycle  is  apparently  one  of  many 
cycles  which  continue  with  a  periodicity  of  roughly  400  ka 
back  into  the  Miocene.  Dissolution  on  this  time  scale 
resembles  the  orbital  eccentricity  record,  which  has  strong 
100  and  400  ka  components.  The  fact  that  the  400  ka  cycle  is 
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not  seen  in  shallow  waCer  oxygen  Isotope  records  suggests  that 
it  is  not  directly  related  to  ice  volume ,  but  rather  to  some 
element  of  the  carbon  cycle  affecting  carbonate 
understaturation  in  the  deep  ocean. 
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1.  Program  DISSOLVE;  This  program  accepts  two  stable  oxygen 
isotope  stratigraphies  and  their  isotope  stage  boundary 
positions,  aligns  the  second  signal  to  the  first  using  the 
isotope  stage  boundaries,  and  subtracts  the  second  curve  from 
the  first.  Interpolation  is  done  with  cubic  splines  (de  Boor, 
1978).  The  residual  curve  is  then  smoothed  and  stored  in  a 
user-specified  file. 

2.  Program  ANALPOW:  This  program  computes  the  power  spectrum  of 
a  stratigraphic  signal.  A  matrix  of  downcore  measurements  and 
depths  are  read  in,  equally  spaced  with  cubic  spline 
interpolation,  windowed  and  Fourier  transformed.  The  raw 
spectrum  is  band  averged  (see  Otnes  and  Enochson,  1972)  and 
converted  to  relative  decibels. 
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Spectral  Effects  of  Time-Depth  Nonlinearities  in  Deep  Sea 
Sediment  Records:  A  Demodulation  Technique  for  Realigning 
Time  and  Depth  Scales 

Abstract  (Chapter  7  of  the  dissertation): 

14 

C  dating  and  %CaCO ^  in  late  Pleistocene  sediments 

suggest  that  deep-sea  sedimentation  rates  vary  cyclically  and 

that  this  cyclicity  is  related  to  climate.  Sedimentation  rate 

variability  leads  to  nonlinearity  in  the  time-depth  mapping 

function.  This  nonlinearity  can  have  profound  effects  on 

signal  spectra,  leading  to  the  development  of  harmonics  and 

intermodulation  tones.  These  distortion  effects  in  the 

spectra  give  a  direct  indication  of  the  degree  of 

nonlinearity,  thereby  providing  a  tool  for  realigning  time  and 

depth  scales.  A  tuning  technique  is  developed  which  assumes  a 

direct  link  between  climate  (as  measured  in  0-18  from 

planktonic  foraminifer  tests)  and  sedimentation  rates.  A 

criterion  of  "spectral  simplicity",  as  quantified  in  the 

varimax  norm,  is  used  to  demodulate  the  input  spectrum. 

Application  of  this  technique  to  an  equatorial  Pacific  piston 

core  (ERDC  84)  found  peak  glacial  sedimentation  rates  to  be 

30%  higher  than  peak  interglacial  rates,  a  figure  in  good 
14 

agreement  with  C-based  estimates  from  the  same  area.  This 
technique  is  compatible  with  other  time-scale  tuning 
techniques  such  as  those  using  orbital  parameters,  and,  in 


combination  with  these  other  techniques,  provides  a  method  for 
fine-tuning  any  late  Pleistocene  record. 


1.  Program  UAVGTEST:  This  program  examines  the  effects  of 
systematic  depth-scale  perturbation  on  signal  spectra.  A 
sinusoidal  signal  is  generated  by  calling  WAVGEN.  The  program 
then  adjusts  depth  values  according  to  signal  amplitude, 
interpolates  with  cubic  splines  for  equal  sample  spacing,  and 
calculates  the  power  spectrum.  A  range  of  degrees  of  signal 
distortion  are  examined  via  a  DO  loop  and  the  results  in  both 
time  and  frequency  domains  are  stored  in  a  user-specified 
file. 

2.  Program  ISOJUST:  This  program  is  essentially  the  same  as 
WAVETEST  except  that  it  uses  measured  isotope  data  as  input 
rather  than  calling  for  a  sinusoidal  signal.  The  input  signal 
is  progressively  distorted  via  a  DO  loop,  and  the  power 
spectrum  is  calculate!  at  each  step.  The  spectra  are  band 
averaged  and  a  varlmax  norm  is  calculated  for  each  spectrum  to 
determine  its  simplicity.  Isotope  stage  boundaries  are  read 
in,  which  are  used  to  quantify  the  signal  distortion. 

3.  Program  VARIMAX:  This  simple  program  calculates  a  varimax 
norm  as  described  by  Wiggins  (1978)  for  quantifying  the 
simplicity  of  a  signal.  The  program  accepts  a  vector  as 
input . 

4.  Subroutine  WAVGEN:  This  algorithm  generates  a  multi-component 
sinusoidal  time  series  of  length  NS.  The  component 
wavelengths  are  specified  within  the  subroutine. 
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III.  Program  Listings 
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PROGRAM  ANALPOW :  COMPUTFS  FOWiR  SPECTRUM  FROM  X-Y  PAIR 
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uPeNiSn IT=1 5 ,NAME= INF  ILF .STATUS* 'OLD ' ) 

:*RIIEI6»40) 

T?  ■r*  aa  A  ft  I  9  t  .1  IT  1  in  T  f  »»•»  A  TT  fT  TMT  fT  r  TT  r  H  1MC1  O  '  1 


FORMATt 

IMS  I'll1  oc 

oniUCRIT-lS, 


WHAT  IS  TEE  OUTPUT  FILF  NAME  ? 
41)  OUTFILE 


NAMF=OUTFILE, STATUS*  'NEW') 

S  THE  OUTPUT  SAMPLE  INTERVAL:  F4.1  ?  ') 


0FEN(UNIT=15.NAMF=0UTFILE.STA 
WH  IT  £  (6 . 12 1 

formatT  what  is  the  OUTPUT  SAM 
READ?!, 13)  XINT 
i0nMATTi4.1 ) 

Input"eata"" 

REAt(l5.*T,rEFTHl(I).XIS0(n 
CUBIC. SPLINE  INTERPCLATICN:  N 


AlHTA.Vi.wf  /  i-  Ci  i  lUi  \  1  /  %  A  sJ  ’  i 

CUBIC  SPLINE  INTERPCLAT ICN :  NEW  LFPTHS  IN  IEPTH2 

npts=(Tdefthi(ns))/xintM 
t£ PTH* 0.0  .  „  „  , 

EFPTE2(1)=DEFTH 
EC  30  1=2 .NPTS 

”^FT^2U?-rFIfcNT 

CALL  INTERPOLATE (XISO, TEPTHl , AIS02, EEPTH2, WORK, NS, NPTS, IBAD.NF.NL) 

NCPNT=NL-NF+1 

ANALYSIS  OF  LATA  , 

kkfcspy*£k5Sc!gfS0SiEK  mu  , 

CALL  COMPLEX(AISC2,CWCRK,NCPNT) 

ZEROMEAN  ANL  LETREND  \IT  NECESSARY  TEE  DATA 

CALL  PROLATE(CWORI,NOPNT.l) 

CQMFUTE  TEE  F0URIE$  TRANSFORM 

aB*»H  HHWH1.BQFW) 


CALL  FlT(Cw6l?RVNOPNf  ,MQPNT ) 

CALL  PCWER(CWCRK,MOFNT) 

CB(I  %CMPf,X(0.2,0.0) 

CALL  CONVOLVE ( CWCRX, CH , WORKOUT, MCPNT, 5 ) 
NPTS=f CPNT+4 
CONyRT  TQ  DB 
PEAK*“1 .0E+38 


DO  50  1-1 . NFTS 

50  ||A|^|X|Fj||gRJAL(WORKOUT 
ELSE  T 

WORKOUT (I )=CPPLX (10. *ALCG10 ( RIAL ( WORKOUT ( I ) )/PE4K ' 

51  cI%l|8!  results 

f'P*i«0PNT/2l  +  l 

gALLRR|ALF|!gT(r^TA, WORKOUT  ,KPTS ) 

AK=(AI/(r,OPNT*XINT)) 

I?UlEAf A  L 1+5 ) .GT.-40. )  GO  TO  90 
IATA7i+4}=-40.0 

90  D0UTU.2  =DATA(I+4)  .  , 

li  f§SMIf-sfl?iJ.9'?T<I,J>,J'1,2) 

STOP 
INI 


va'(  1.00  ^.icrnii4'>-NO'i-*cs;u'jnr -orr>oit<s.r^oii-.(^.(oo 


HRiMSoraiflfcTJtlci^"  r*CCNV0Lim0N  ««  rEC0 

DIMENSION  XCUT (166) 

CfaAKACT£R*2e  INF ILE1 .  INF  IU3 
CHAHAC1ER*20  CUTFIlFl 

*RITE  (6»1 

1  FORMAT  .NUMBFE  OF  SAMPLES  IN  INPUT  ISCTCPE  STRINGY:! 
1»RITE(6,3) 

3  FORMAT  l  * .WE AT  IS  TFE  ISCTCPE  INF I IE?  ') 

4  Httiitl)  INFIlil 

CPFN  (CNIT*1S,NAMF=INFIL51 ,ST»TUS=  'CLD') 

MP1S=64 

...  iwa.i-1:64 

11  FCRKItT VVHAT  IS  TKE  ISCTCPF  CtTFIIEV  ) 

RFAI  (5,4)  00TFIIF1 

Of  SN  JyNITj^e.NAM-CUTIHF:  ,ST  ATUS=  NFW  ) 

12  FCRMATi*  ¥  :■  A  T  IS  TEi  UCC  n-  >■  I  \  LIMIT?  :F17.S  ') 

READ  £,13)  G  *  I N 

13  FORMAT  til  </  .9) 


52  mmusuhih™  ut  !':n!i " 

C?EN(lMT=25,NAM  =  lMlLE3.STATUS  =  'Oil  ') 

r. i AI  IN  LATA  *' 

DO  22  1=1 .  NFTS 
22  SEALTiS,*)  XISC(l) 

NE  =  NPIS-*r  ETS-1 

CALL  rECC(h.NPTS,XISC,vClT ,04 IN 
10  34.1-1. Nf  ,  v 

ii  mmSyomi) 

STOP 

INI 


19 


12 

22 

22 

42 


72 

82 

se 

120 

112 

120 

120 

m 

170 

180 

192 

222 

212 

220 

222 

242 

2t2 

262 

270 

282 

292 

222 

212 

222 
T  1  Q 

Jy  V  1C 

242 

250 

262 

270 

280 

2S0 

400 

410 

420 

420 

440 

450 

462 

470 

482 

(if 

522 

522 

ill 

560 

572 

^82 

!se 


C 

C 

C 


f h8g!ak  a{3istseiwo  stratigraphies  about  input  control  tips 

ANE  CALCULATES  THE  RESIDUAL. 

DIMENSION  XXXI  (200) .XXX2  (200  ' .  XXDl  (200)  ,XXD2(200  ) 

(30),  FACTOR  (30)  tHjl07 

Stall  8K  . 

\FILE1,}n!1L£2 


CHARAC 


200)’XISO2(200),XOUT(200) 


ers  im: 


CHARACTERS  OUTPILE 


1 

2 

3 

4 

7 
99 

8 
c 

10 

11 

12 

13 


WRITE (6,1 ) 

FORMAfT "FILENAME  FOR  FIRST  CORE? 

iikw,i‘‘ 


CERAllf J 

7.1=1 


10 
RE 
CA 
FORMAT 


14 , NAME=1NFIIE1 ,STATUS=  'CLP') 

CW  MANY  SAMPLES  IN  THE  EIRST  DATA  SEQUENCE 
NS1 


:  I  2  ? 


XXDl  ( I ) .XXXI ( I ) 

LL  ZERO  (^ XXI, NS  1. /.ME AN1  ; 
fvRMAT(5X,F5.0,3X,#5.27 
WRITE (6, £J 

FORMAT  (  \E0'*  MANY  STACE  BOUNDARIES 


NSTG 


REAr(6.9) 

FORMAT ( 12; , 
WRITI(6.10) 

FORMAT ( *  ' 

DO  11  1  =  1 _  % 

RFADT5  ,*)  STGl ( 1  ) 


? :  12  ') 


STAGE  BOUM ARIES  IN  CM  ') 


JoRMjHf  ^LFNAMF  FOF  TEE  SECOND  CCRF?  ) 

RFAI (5,2)  INFILE2 

0FENiUnIT=15,NAMF=INFILE2,STATUS=  'cir  ) 

WRITE (6. 13) 

FORM AT (  HOW  MANY  SAMFIES  IN  THE  DATA  SEQUENCE?  : 13  ') 

D0*1*  I  i  NS" 

FORMAtTF4.0,F6!2) 

WR  ITFlc  *16) 

FORMAT i  INPUT  STAGE  BCUMAPIFS  IN  CM  ' 


71 

73 

72 


30 


ANr  MULTIPLICATIVE  FACTORS  ') 

mnvz:mi{iVs™2h) 

FACTOR jl) = (STGl (I)-STGl (1-1 ) 5/ ( STG2 ( I ) -STG2 ( 1-1 ) ) 

MtifeSItfS1 » 1  •  sic?  ( !  ■ .  f>.ctob:  i ; _ 

I!i&)%?FfiETHE  cu,put  mt  ■'  ■> 


•J 


OFEN  (UNIT  =  ie,NAM=OUTFIIE  .STATUS  - 'NEW' ) 


STANDARDIZE  SAMPLE  SPACING  IN  CORE#l:  10CM  I NTERVAIS 
GENERATE  A  NEW  IEPTH  SCALF 
NPQINT=XXD1  NSlT/10. 

4s  WHlMoR?!4*.!*) 

xc(i]=e.0 

25  xS(i7=x5l:i-i?li0.0 

CALL  INTERPOL  A.  TE1XXX1.XXL1 .X  ISO  1 .  XT . WORK .NS  1 , NF01N T, I  BAD, NF, 
REFERENCE  ISOTOPES  ARi  IN  fal.XlSoi);  StfiOM  SFCUFNt'E  IS  IN 
(XXD2.XXX2).  NEW  DEPTH  VALUES  FOR  CCRE02  ARE  FIT  IN  XNEW. 


TO  110 


J=1  ,  , 

XNFV(1  )=0 . 0 
10  100  1=2. NS2 


WE  ITII6 .42  I ,  XXE2 ( I ) , . 

IE  IXXEJTi  .l£.STG2( J) )  GO 


WRlTE(e»432 
43  FORMAT (  NEXT  J  ) 

1  le  1 XS  If  1 ?*Di  !  1  I‘>  >  )*F»CTCS  U  ' 

42  ICRMAT(iS ,iexlFl0.i; 

100  CONTINUE 

INTFRPOLA  TF  T  V  A  LUES  FOF  NFW  PIFTFS 

CALI  INTERPOLATE (XXX2.XNF/ t  XIS02 . XI . WORK , A  S2 , NFG IN  T , I B#  D, NF, 
CALCULATE  RESIIUAL  AN f  PI  A.Ce  IN  $RE §  IT  . 

DO  200  1  =  1 , NFCINT 

2ee  imvrnHimhxmn}) 

MPTS=5 

F(1>0. 11111 

r  \c ,  ~a.cc.ccc 

BC3J=e.  33333 
6(4] =0.22222 
E(fJ=0. 11111 

CALL  CNVLVjXRFSII .K.XOUl J POINT  .NETS 
NPTS=NFCINT+MFTS-I 


10  201  1  =  1 , NPT S  ,  t 
WRITE  (10,  202  *  XDd).XCUT(n 
FORMAT (E10.1.10X ,Fl?.  5^ 


STOF 

EM 


I 


PROGRAM  GSTEST  (GUINASSC-SCEINK  TESTER ) 


LIMENS ICK  A (64) 

CHARACTER*??  OC1IIIE 

i £ j 

10  IDEM  AT v  *  WHAT  IS  THE  OUTFILE?  ') 

RIAE(5,12J  CUTFIIE 

1 1  FORM  AT  ( A  ' 

open  Tun  it  =i  e  t  n  ame=outf  he  .  sta  tus=  'new  ' ) 
WRiTElene; 

20  FORMAT!  WHAT  IS  TEE  MIXED  LAYER  IF ICK  NESS  :F4 . 1 
REAL  (£  ,21 }  LEE 

21  FORM  a!  If  4 . 1 ) 

WRITEt6*E2T 

22  FORMAT!  WHAT  IS  THE  riRETSICN  COEFFICIENT?  ) 
READ  (£.2l!  Lit 


WRITEjS*?^) 


AT  IS  TFE  SEE  PATF:F4.1  ? 


31  FORMAT  (14.1 

CALL  GSIRFTA.IEP.CIF.V 
DO  101  1=1.  el  , 

WRITE! 16,1  A  I 
101  WRITE  6,1  f  ( I  v 
1  iORMAT  U0A.I6.6' 

STOP 
EM 


11  FORM 


f’^gglER  OF  SAMFXES  IN  INPUT  ISOTOPE  STRING?:  13  ') 
'$! 

l*.WEAI.IS  THE  ISOTOPE  INIIIE?  ') 


.41  INFILI1 

NIT=15,NAME=INFILE1 .STATUS* 'OLD ' ) 
?*^WHj}JiJ|2THE  ABUNDANCE  INIIIE?  ') 
AlT*25, NAME*INFILI2 .STATUS*  'OLD ' ) 


4 I* 1.64 
•  ® 


f*^AJFIJ  THE  ISCTCPE  CUTFILI?  ') 

|A!t*?8, n1k1=CUTI IIE1 .STATUS*  'NEW  ' ) 

I^MHAT  IS  THE  EECCN  GAIN  LIMIT?  :I10.9  '> 


nPONSE  FUNCTION 

AJiI|3THF  IRI  INIIIE?  ') 

5, NAME*INFILE3 .STATUS*  'OLE ' ) 

51 

DO  20  1*1 .NPTS  ,  . 

« w*i 

ha  M;s«i  inriiw’ 

2?  aintia^MD 

26  MIf|<f*f^pA4(1) 

SMita) 60 10  46 

is  lu¬ 

ll  IB«m8i?®i)*5(I> 


£6  Si m (I 


GO  TO  40 


A5(  I ) 


( 


ON  fcITH  SIGNAL  CARRIER  ABUNLANCE 

|i.E{647 

A3Tl50j.A4il50)  ,A5  150) 

IE2  . INI i LEo 


15C 


i/iiiLiiu  ivn  nil  /  •  nv.  V  1../V 

CEARACTERsi‘20  INFIIE1.IN 
CHAR4CTER*2e  OUTFILEi 


1  I0RMAT(  NUMBER  OF  SAMPLES  IN  INPUT  ISOTOPE  STRING?  :12 
READ(5«2)  NPTS 

2  FORMAT? 12) 

WRITE(6,3 

2  FORMA! ( *  WEAT  IS  THE  ISOTOPE  INEIIF? 

REAI  (5,4 )  INFILE  1 
4  frORMAfUj 

CPE  N?UM1  =  15,N  AM.  F  =  INFIIE1. STATUS®  'Oil 


i:RITM6,5) 

5 1?SE1iT4)v?£i:iIi-THF  ABLN1ANCF  IKFIlF? 

OPEN  (UN  IT  =25  fNAM,E  =  INEIIE2,ST£TUS='OLI ' ) 

MfTS=64 

DC  102  1=1.64 

100  nTi]=0.e  „ 

WE  I TF (6  ,1 1 ) 

11  FORMAT  ;\»cAl  IS  1  Hi  ISCTC?.1  CU1FIII?  '} 
READ  6.4)  OUTFILU 

CPEN(UnIT®16,A  AMF=CUTF llFi .STATUS*  NF*  ) 


IMFULSE  RESPONSE  FUNCTION 
WRITf?6,50j 

52  FORMAT  (  WH AJ  IS  TIE  IFF  INfILF?  ) 
EEAD(*,4)  IN/IIE? 

OPIN  (UN  IT  *35, NAME* IN FILES .ST  ®TUS=  CLL  ) 
DO  51  1  =  1  *  64 
51  E  E  A I  (  5  5  ,  *  J  H(I) 
fi  EAit  IN  LATA 
IC  22  I  =  ! , NPT S  , 

22  READ,U5.*  XlSOm 

10  2<  1=1, NPTS  , 

22  REAI  (25,*)  EOEd 
DO  ,2  5  1*1. NFT? 

25  HUM U  pcpn! 

CALL  CNVLV ( FOr  ,F  . A4 . NPTS . MpT S ) 

CALL  CNVLV  A1  ,fi  ,  A2, ms  .MTS  ' 

LO  34  1*1, NP 

11,(14(1). EC. 0,0)  GO  TC  4/ 

A5(I)»A2(1 )/A4(I 
GO. TO  33 

3§  cSAtinEe0 

34  WRITEjlc.26)  A5H 

26  FORMAT ?Gi7.r) 

f&F 


24 


I  i if  in  n(  T«c  w  <  i  I  f  III  m  no.-.  a-.  **r*iC  *fji  Mf-. »i.',.r.,r  J.  •.Mn-lvi'*  'MMn  ■>-  *■  •  • — »►-— . — >  > — » »— • . 

uim-oOii  Ha»f»ifot-»MiiiiJ  -»xt>i  maruM^jiod  -oi">-i  utp  -h  hi iiu.>'mai»-'<m(i  il  ->xm  -vicr.il: u»i  -imm 

cviSv <S:r>  fvw:rjfv;  ,{S;e:;f^v^s‘.<^P-e:;MPsJ<t:<'v  e  _■  f'.is.  c_  i^.f -ini  <v-  mvis  r-  ^.is.-r>ocsj  i\vP-.;csi-.  <r:  rcivnr  in.  c.  inin/'-' ivin  inis; 


c 

g 

c 

c 

c 


PROGRAM  ISO JUST  (ISOTOPE  ATJUST) 


PROGRAM  RECALCULATES  ISOTOPES  WITH  DEFTH/TIME  ADJUSTMENTS  BASED 
ON  S  FDIMENTATICN  RATE  CPANGES .  DEPTH  SCALE  IS  MODIFIED  LINEARLY 
ACCORDING  TO  SIGNAL  MAGNITUDE. 


DIMENSION 

DIMENSION 


£Im5NSI0N  XISQ(50 

DIMENSION  XN3T500 _  __  „ 

COMPLEX  CWORX (500) ,  WORKOUT (f 
CHARACTER*20  INFItF 
CHARACTFR*20  OUTFILE 


.DATA 
2 ] , COU 
Hi (500 
00), ST 
CK(10) 


(500,2) 


60 

1000 

51 

40 

41 

125 

126 

12? 

12c 


WRITE (6 *80) 

FORMAT  1  INPUT  NUMEER  OF  SAMPIES  IN  THE  DATA  SF0UENCE:I3  ) 
READ  (5.1000)  NS 

,0SE 

FCRm|t(  yHAT_IS^THE 


.  ..  T_.  v  ______  INPUT  FILE  NAME? 

RE A 1(5, 41)  JNF ILF 

OPEN (UN IT  =  15. NAME'INFIIF, STATUS'  OLD  ) 
WR IT  m6  .40) 

r  r\  n  »  m  ■  ~  i.i  rt  1  in  T  n  m  r:  r  r 


THE  OUTPUT  FILE  NAME  ? 


) 


FORMAT  (  -  _ _ 

till5  NIS0 

HP.  ITE(6,127) 
FORMAT l  CORE  DE 
READ (5,128)  ISTG 
FORMAT^ 20F4.0) 


OUTFILE, STATUS'  'NEW') 

OF  ISCTCPE  STAGES  IN  SERIES?  :  13 


FTPS  CF  IS< 
v I ) , i=l  ,ni: 


'T^FE  STAGE  BOUNDARIES?  :20F4.0 


') 


_  c 


80 

50 

501 

56 

81 


READ  DATA  INTO  XXX 
DO  52  1=1, NS 

lit;  M)'-2’ 

NFWl!l|prS(H|h0,n' 

ECRMAir  NPNTS *  ,14 

UATE^DATA  INTO  10CM  INTERVALS 


INTERPOLATE  DATA  lUi 

CALL  RAN'GETxISO.NS, 

WR ITE(6  *501 ) ,EX  ' 

DEPTH2(l]*Il-l 
CALL  INTERPOLA 
NPTS®NL“NI+1 
WRITI(6*91)  NP 

Tnnti  iip  1  *  iutfc  = 


DMAX.DMIN ,DX ) 


‘U 

E(XIS0,DFPTF1 ,XISO 2, DEPTH 2, C WORK, NS, NPNTS, IBAD,NF,NL  L 
n£tS'T5,I4) 


FORMAT! 

ZFRCMEAN  DATA 

EltTLH2AE!3ilfcNPTS’^UN) 


25 


J  t  t. 

zZL 
J  C 

ySU 

\<ll 

Ul't 

li'il 
x  c-i  £ 
\i<kl 
1  /52 

\  xZ  < 

li c*U 
12=0 
Vil'i 
il  It 
lilt 
1  lk  2 

,  i.7/ 


12  =  0 
li'?'7 
1  1  £% 
1 1 9  v? 


ANEWU  1=0.2 
I)C  73  1=2  NPTS 

xNEW(I)=C0UNT+(ie.6>(l.-TER^) 


GC  TC  73 
xnewTi  =C0UNT 
CCUNT=ANLW(I) 


+(ie.0)*(l.+TERN ) 


WRITE 

FORMA 

1=1 

DC  112 


ETERNINE  EFFECTIVE  SEDIMENTATION  DIFFERENTIALS  ') 


DC  1 12  1=1. NPTS  ,  v  x 


L  =  L  1 

IF  (I.GT.MSO)  GC  TO  750 
CONTINUE 
CONTINUE 

a' P I T E  CUT  AEJl’STFD  ISCTCPF  STAC-E  BOUNDARIES 

FORM »Tr  VSjUSTFD  ISOTOPE  STAGE  BOUNDARIES  ' ) 
TO  113  1=1, NISO 

EFFECTIVE  *A 


WPITF(6j114)xIjXST<^(I  ) 

EFFECTIVE''*  ACTCR  *  ( LENGTH  OF  ADJUSTED  STAGE )/( LENGTH  CF  INPUT  STAGE) 
•vti  IT  E(6,lie) 

FORMAT  t  EFFECTIVE  S  ED  IMFN'T  AT  IC-N  RATE  FACTORS  ') 

M|cT=(  (XSTG(  I))/(STG(I))) 

*?  II F  (c .  1 17 j  I.EFACT 

FORMAT  [  ISOTOPE  STAGE  =  M3,  FACTOR  =  \F5.2 
L'C  118  f  -2  NT  sc 

,i:ms4I-1,:/(STC(i:-ST5<1-i),) 

CONTINUE 

SJCEEZl^OR^SJUTC^y^EW  TO  SAME  LENGTH  AS  XIS02.  (  =  NPTS) 

WRIT  tie *10)  NF 

FCRMATr  NUMBER  CF  POINTS  =  ,14) 

ANPTS=NPTS 

ANP=NP 

^IHiir»H8iMi,EiN8fF"SiTitt!UTi  yITH  THE  m  IEPTH  SCUE 

,  RECALCULATE!  ISOTOPES  ARE  CONTAINED  IN  XN4. 

LEPT  R=« .0 

IHJ’iU  SPts 

DEPTH  =DEPTH+10.0 
XN2(I)=rEFTH 

CALt  INTEF|01ATE{XIS02,XN1 , XN4 , XN3 .WORK , NPTS , NPTS , IBAD , NF , NL ) 
NOPNT=N  L-NF+  1 
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1222 

:zu 

1222 

i  Z  1 
Yell 
1262 
1272 
126  2 
12  2  c 

I  ill 

:  liz 
1222 
1232 
1242 
lie  l 
12c  2 
i  ?  ?  2 
1262 
12s  2 
1422 
1  112 
1422 
14  52 

I I  4  o 
1452 
142  2 

1 4  f  2 
1462 
1442 
U22 
1  c  1 2 
1522 
1222 
1242 

1 5  5  '2 
166  2 
1272 
i5c2 
1552 
1522 

im 

1621 

1622 

1620 

lile 

1672 
1662 
1692 
1722 
1712 
.  1722 

*  1722 
1742 
1752 
1762 

•  127? 


MBlPoF  POINTS  ARE  =  ',13) 


CO  102  I  =1 , NOFNT 

Wjttfem,,,.,, 

FORMAT v  1X.2G15.7 ; 

CONTINUE 

ANAiySlS  OF  DATA  .  _ 

CALL  CLEArt CWORK.500 ) 

TRANSFER  LATA  INTO  WORK  AREA. 

misimmwhyhihs>. 


mirm  fETftERr  TrrsrtfssAsi)  thf,  tm 

CALL  ZEROMEANTCWCRK.NCPNT) 

*  I  MO  Vi  THE  TATA 

CALL  PROL AT EfCWORK  .NOFNT .0) 

COMPUTE  ThE  FOURIER  TRANSFORM 

mopntmbinaryTnopnt) 

CALL  FFT(CWORK  .NOPNT  ,MQPNT  V 
CALL  POrfEPlCWOfiK »MCPNT ; 

PANE . AVERAGING  e 
Cfi  '.I  )  =  C.VPLA(0.25,0.  ) 

CHC2  j  =  Cf'PLX(0.5i,2.  ) 

CAU  CCNVCnflCW{RK!cR.'*CRKCUT.MCPNT,3  ) 

CONVERT  TO  RELATIVE  EEC  I  BtLS 
PEAK=-1.0E*38 

192  f  FAMMA|7HS^j|AL(IVCRFCUTf  I  ))) 

IF  [pEAL( WORKOUT { I ) ) .EC .2 . )THEN 
^R|gu||l)*Cr'FLx(-l00.,2.) 

EL3o*KOUT(I)=CMFLX<10.*ALCG12(PEAL(WORKCUT( I )  )  /PEAK)  ,0. ) 
EM  IF 

191  CONTINUE 

OUTPUT  RESULTS 

CALL  SeaBaRT(DAT.A  .WORKOUT, MOFNT+2) 

AIM-15  *s<:,r*  +  1 

AK=AI/ (MCPNT*10. )_  __  _ _ _ 


AT  A ( 1  + 


45  l*RITE( 
47  FORMAT 

VARIMA 

&IJ 

V=2 . 2 
DO  95 
V1=V1+ 
V2-V2+ 
95  CONTIN 
V*Vl/( 
'•RITEj 

46  FORMAT 


T .-25  .  )  GO  TO  45 


wimf jltii)  AK*fIATA(I+l) 


3??i;sWw.- 


27 


i?se 
it  ze 
IS  le 
la2fc 
le2e 
1 E  4  2 
late 
ItcZ 
it7e 
id5e 

lesZ 

is  2  e 

ISlz 

Is  2  2 

i53e 

1S42 


CONTI 

S88I1 

ENE 


yi8£y2I$!?IrR£N5£lJC*NS.Xf'AX,XMIN1DX) 

IMEnI ION ElF§Sl¥^S  rA**,'lK  VAIU^S  ANE  AMPLITUDE  OF  A  VICTOR 


DIMENSION  U500 
XM IN*X ( 1 ) 
XMAX=XMIN 

?9  f X?I ?7lf ^Imi N 

ez  HH  HiGI-xru 

IX=XMAX-XMIN 
RETURN 
ENt 


XMIN»X ( I ; 
XMAX=X ( 1 1 


ffHi 

.THE 


5o58itF5! 


iU  uniniil  lEi  rii  V,r 

IPS  P M«W»itnBiwiiWw> 

DIMENSION  CUT3(I28) 

98l8M:HiB»ffltiii,u,,-l"lt,8-1"11M,1MI1K _ 

up itt ( g  220) 

FORMAT v  HOW  PANT  I ATA  POINTS  :I2  ?  ') 

nlmmy Np 

BEAD  IN  DATA  FILES 

BtflMil  FILENAME  ?  ') 

READ(*203)  INFILE1 

£?EN (0NIT=15,NAME=INFILE1 ,STATUS= 'OLD') 

IC  204  1=1. NP  ,  . 

CALL  zi§C?AlSCl^N?JpFAM  ' 

Ais§??il=ij,soi(i)*0.8? 


FORMA 

WRITE 


m 

.206 

W  T  M  A 


OPEN  (UrtlT=4?,NAME=INE  1LE2, 


STATUS-  'cir' ) 


?  ') 


ur  in  i  uiu  i  -  i  r  •  w  ji  •  •  r  ”  a  w  *  a  i-  a-  *•  «  v  .• 

Ca{i  flfto (XlSOS.Nt .APE AN?' 

FCRMAtI '*pf-Pl  FILFNAPE  ?  ') 

c|en!p!t*18.  MPF  =  INFUE3  .STATUS*  '0LL  ' ) 
DO  209  1=1. NP 
READ(ie,20fe)  POPl(I) 

ifcRMAU  'fifyz  rilEKAKF  ?  ') 

oIen  1ui5iT=19!nAPE=INFIIE4  .status*  'CLI ' ) 
ro  211  1=1. NP _ 


f?ASU9!|JfTF0F2(I) 

I0R«a|! ^WEIGHT  FILENAME  ?  '} 
olEN(pIT=l0.NAp|=tNFILE5,STA 


ST*T0S*'01B') 


9AlS*l*1J  iFSIDlATI 


ION  CHANGES 


,a  iMill 


I 


29 


DII-  '.F10.5.5X,'  E  =  '  ,F1k5. 


read  in  iaiTIAL  VAICES 

FORMAT v  >'HAT  ARE  THE  INITIAL  VALUIS:EEP,IIF:2F5.P?  ') 

bsImW’  ■ 

FCRMATl*  what  IS  TEE  SFI  RATE  :F4. 1  ?  ' 

1SMI  pM 

WRITi(6.13r 

format!*  WHAT  IS  TEE  DECON  GAIN  LIMITATI OK :F10 .9  ?  ') 

ipplSV" 

CALL  ERROR  X,E,V  ,PU  1  ,PU  2  ,P*1 .  PW2,  GA I N.  NF .  OUT1 , 0UT3 
WRITF(6,12)  i  ( J ) .X (2 j , F  ,  , 

FORMAT ( *  DEP=  ,F10.5,5X.  DIF-  '.F10.5.5X,'  E=  '.Flrf. 

INITIALIZE  SIMPLEX  ~~  ~  ~~~~~ 

DO  22  J=1.2 

fN’!& 

lilAiiHri-lTi.IiJiE-iE,)  ? (i,i-i )*0.000i 

CRIER  POINTS  ACCORDING  lO  V*IOES  CF  FRRCR 
L-l 

ft  M-3 

CAU  FRROR(X,F.V.FVI  l.PkIP.FVl.  PWP,  GAIN,  NP  ,  OUT  1  .CUT3 
fluRslf  ).LT.ERR(L))  L  =  I 

if(err(i).gt.err(hM  K*l 

NE=L 

LC  43  1=1 .3 

IF  (ERR  (I)  .GE  .ERR(NH  )  .AM.I.NE.E)  KH*I 
CALCULATE  CENTROID  (AVERAGE) 


Mar 


FORMAT J^RjpLECT  AEOUT  THE  CENTROID') 

Cm*  !v?f8??  Jpwlljpil  ,PW2,GAIN,N’P,0UT1  ,CUT3  ) 

!!  ,a8°tt*tlnr  sctcISSFOl) 

IF(ER.GE.IRR(H))  GO  TO  122 
DO  60VJ=1.2 

|I|.J  -Rp) 

IHER.GT.ERR(NH)  )  GO  TO  61 


.ERR(NH) )  GO  TO  61 


12  le 

liiS 

1240 
1250 
1 2ee 
i2?e 
1280 
1292 
1222 
1212 

i  m 

1242 

i  m 

1370 
1380 
1292 
1422 
1412 
1 421 
1430 
1442 
1450 
1462 
1470 
1480 
1490 
1520 

1111 
1122 
H42 

1112 
1?60 
1670 
1680 
1590 
1^00 
1610 
1620 
1630 
1640 
1660 
1660 
1670 
1680 
1690 
1700 
1730 
17  20 
1730 
1740 

i77e 


WR?TE(6*503 
FORMAlI*  IX 
E=NE 


IXl  AM  THE  S  IMPLEX  ') 


DC  93  J=1^2 

II (EX. it. ER)  GC  fo  104 


imi* 

GO  TO  110 
DO  10*  J=1 ,2 

pTi,  j)=x( j 
err{l)=ex  % 

WRITE  6,12) , 


12)  P  (  L ,  1 ) ,P (1,2) . ERE ( I ) 


WRITI(e;ill) 

RIAll*  11)  I?¥  *ISE  TC  CCN1INDE?  YES  =  1,NO=0  ' 


READ  (5 , 112 )  IEI 
FORMAT (II)  % 

II  IIE.E5.1  GC  TC  41 
GO  TO  150 


WRITE (6 *604) 

FORMAT  (*  CONTRACT  THE  S1MPIFX  ') 

R  ( j)  =2 . 50  l£*C  ( J }  ->  0  .4985 *P  •'  E ,  J  ) 

CALL  ERROR (R.iR.V,P*Il.EWI2,PVl,PV2,GAIN,NF,0UTl .OUT 

n  |»!Ei‘Ki.h1  Msfii 

II  (EH. LI. ERR (h))  GO  TO  $9 

SCALE  TEE  SIMFlEX:  SK  IS  TEE  SCALING  FUNCTION 
WRITE  (6»130) 

FOPMATl  WHAT  IS  THE  SIMFTEX  SCALING  FACTOR  :F4.1  ? 
READ  15.131 )  SK 
FORMAT  v 14 . 1 1 

p?i!j%pu!j)+sk*(p(i  ,  J)-P(I  ,  J)  ) 

GO  TO  31 


TRACT  THE  S1MPLFX 


!j)+SJC*(P(I,J)-P(I,J)) 


iriTF^tlSl) 

FORMAT r  WHAT  IS  TFF  CCTIIIE  7 
READ {5,203 )  OUTFILE 

OPEN  (UNIT=16,NAM=0UTI  IIF.STATUS=  NEW  ) 

NPTS=NP+64-l 

A= CNFTS/2 J-20 

B=  (NPTS/2  j+20 

DO  1 52  I=A , B .  4 

ouiiTn*ouf im*i.i5  % 

WRITE(16. 153)  I. CUT1 (I  +AMEAN 1 .0UT3 { I  +AMEAN2 
WRITE (6. 1531  I  .CUTlfl  ) +AMEAN 1 ,0UT3( I J+AMEAN2 

iormatuS,  iex,di?.5,i0x,Gi?.5) 

INI 


t 


c  PRCGRAK  PROJECT:  CCI* 
DIMENSION  X 1 K Tl2 


CTER’ 


W EP  .SPECTRUM  IK  L'B 


NFlIF.Ot'TFILF 


fcR  IT  E  (6  1 ) 

1  |^5‘J^2)NEffsR  °F  SAf'PIFS  IN  INPUT  :I2  7  ^ 

2  FORMAT tl2)  F 
f«R  IT  I  (c  3  ) 

3  siSM&i'mu  tee  iniue?  ’> 

4  FORMAT Ik) 

OFFN  jDK' I J=15.NAME=INFIIF. STATUS* 'Oil 
£  FORMAT? *5WFAT  IS  TEE  CUTFIIEV  ' ) 


THE  INEILE?  ') 


£  FOPM 


WFAT  IS  THi  CLTFILEV 


OPEN  (gf5lT=ie,NAMF,=CUTFlLF.  STATUS*  NEW 
1*R  IT  eTc  ,6 ' 

e  FORMAT r  WHAT  IS  TFE  SEE  RAT1:F5.2  ?  ) 

HtfiiiJi.I) 


BEAD  IN  DATA 

10  REAi^i|~il1PxfK(n 

11  F0RMATji£A,F6.5.‘ 

GAEL  CLEAR JVCBK.126) 

CALL  CQMPLmXIN.WCRK.12S 
CALL  FFTXWCRs.Kpis.lSr'F. 
CALL  FOteERDB(WCEK,i?b] 

rc  izz  1=2.65 
x  i  =  i  - 1 

FREC=X I/C 1 M 

U2  V.P  ITiT  16 . 122)  FREC  .XCyT  1 1  . 
122  FORMAT (F1Z. 5. 10X .FlZ.c' 
STOF 
ENL 


7  fct£ 

eze 

see 

3SSS 

1111 

i4Ze 

i*ee 

leze 

i?ez 

18ZZ 

l-zz 

2ZZZ 

2122 

222fc 

22ee 

24  2  2 
2522 

2eze 


*NHlMMMHWIHmcm  K0E” 

CHARACTERS  TNI  IlE 

Ifl5f^2)EPW  MANY  SAMPIiS  IN  TRI  INPUT 

lORMA? fl2) 


2  MMiifK) 

^  m&M 


2  FORKAll'TUHAT  is  the  ineiie  ?  ') 

READIc .4)  INFIII 

OPEN T^Nlt=15,NAME=INFIIE  ,STATDS=  'CLE') 

S  REAlHlsM  XIN(I) 

V  1=0  .id 
V2=0 . 0 
V=2.2 

EC  11  1=1, N  x 
V1=V1+(XIMI)**4) 

11  V2=V2+(XIN(I ) 

V=V1/(V2**2) 

MITEje,10&)  V 

ize  format!  tee  variwx  norm  IS  '.Gl£.?) 

STOP 
ENE 


r  w  tx  ntm  i»<  w  mm  in  n>&  a*  it.it.  rf»itntntir.ic  jcjcjc  .if  jc.jr.iriK  jtiin'tviMJsiM  mmmmmhmm  •-•»-»  ►-»  t-**-1 1-* 

uvn-ofn*  ni^f  jNM-»^s<na.'-ciow  mt»c>ir\ii-‘»«otoG.  ->in\  ( iiii.gjm  >-»  ffitorr  -oou  mt»  r  jm^ps  tooc-vinw  ihiwr  jMHNmm-om  nu^  c  jrvt-1 

F3^3FJSKiF»WK>»^53«^)FCFSK>evr>^««CS>l'.K5afoeN;t,5«,0<l5«^)"3K)ffi(t-K>^  (<OPi  Pv««aN!<^PSea<,3<^><S««3<SWC5<?ae<5C5 


PROGRAM  VAVETEST 


PROGRAM  GENERATES  A  MULT I-WAVELENGTH  INPUT  SIGNAL,  AND 
SIMULATES  SEDIMENTATION  RATE  CHANGES  WITH  AN  AUTO-MODULATION 
BASED  ON  SIGNAL  AMPLITUIE .  THIS  FROGRAM  IS  SIMILAR  TO  PROGRAM 
iN^tll’  1lH*CE  ACCipls  CXYGEN  ISOTOFE  DATA  AS  THE  CONTROLLING 

DIMENSION  XXX(500.2).F( *00j~iuF7|00 7 ~~7rATA 7 5007 

COMPLEX  CWORK  500] 

CHARACTERS  OUTFILE 


1000 


WRITE(6.£0) 

FORMAT!  INPUT  NUMEER  OF  SAMPLES  IN  THE  DATA  SEQUENCE  :I 

HMiiiHPV* 

*R  ITE(6  .40 ) 

T  A  T>  4j  A  IT.  I  ’  T  '  T7  •  m  TO  rr  T"  n  ^  tt  fr  r\  r,  rr  ^  v  T  r-  i  .  i  >  a  /  V 


iter 


ATr  WHAT  IS  THE  OUTPUT  FILE  NAME  ?  ') 
(5,41)  OUTFILE 


NAMF=CUTE ILE, STATUS*  'NEW  ' ' 


GENERATE  RAW  DATA  , 

CALL  WAVGEm  XI  SC  ,NS  ) 

CAIL  RANGE  (X ISO. NS  ,IMAX  .  LM.IN  ,  TX  ) 

CALL  MEAnTxISC.NS.DMFAN) 

WRITE(e,l5)  IMAX  .rMIN.rx.lMEAN 
;  FORMAT!  MAX, MIN. RANG* ,  AND  MEAN 
14 (2A.F8.2 ) ) 

ioAr  xisc  into  xxx 

DO  50  1=1. NS 
XXXTl.lMl-l)*ie.0 

i  xxx(!  ,2Miso(n  ,  , 

NPTS=(((XAX?NS.l )]/10  +1 
WRlTi]6,90)  fcFT$ 

I  FORMAT!  NFTS-  .14 
DEPTH  ADJUSTMENT 
DO  73  J=1 ,6 

W4(i‘1,/ie-e) 

XNEW(1)=0.0 

UeWi]=COUKT+TI0.0)* ( l.-TERM) 

CCUNT=XNEW(I) 

GO  TO  72 

CONTINUE 

SQUEEZE  0£  ST^ETCF  XN^W  TO  SAME  LENGTH  AS  XXX.  (=NPTS) 

WRITE (6 *10)^  '  '  +1 

FORMAT T*  NUMBER  CF  FOINTS  =  ',14) 

ANFTS=NPTS 

icpiiPi=i,NS 


VALUE  OF  ISCTCFE  CURVE 


m 

eze 

tiz 

y  i 

672 

m 

722 
712 
722 
732 
74e 
75* 
7  ee 
7*2 

7e2 

7S2 

822 

e;? 

be.  2 

632 

842 

§1* 
t?e 
86  2 

Iff 

s  12 
g22 

s32 

r'-.y 

see 

£62 

£72 

562 

1111 

fliS 

1232 

1242 

m 

1272 

icae 

1292 

1122 

1112 

HIS 

1142 

1152 

nil 

1170 


c 

c 

c 

c 


11  XN 


EE 


XNltn  AND  Xiscll)  ARE  USE 


C 

C 

C 

C 

C 


E 

_ .  .  .  ...  USET  AS  INPUT 

XN3,  01  LENGTH  APIS  . 
ISOTOPES  ARE  CONTAINED  IN 


WITH  THE  NEW 
XN4. 


30 


62 


120 

121 


M3) 


SPECIFIED  AS 

upturn 

XN?Il  =EEFTE 
DO  30  1=2 ,NFTS 

iKIlirBItti*1*-* 

call  Interpolate  x iso, xm,xn4,xn3, work, ns, npts.ibal 

N0PNT=NL-NF+1 

E§E?A}f >6SoMbSrNJp  FOINTS  APE  = 

STORE  TATA 

EO  102  I  =  1 . NOFNT  ,  v 
WRITE  16,101)  XN4  I 
EORNATflXlGlS.?] 

AN  A I YS IS  6F  DATA 
CALL  CLEAR (CWORK, 500 ) 

TRANSFER  DATA  INTO  WORK  AREA, 

CALL  COMFLEX ( XN4 ,C WORK  .NOFNT ^ 

ZEPOMEAN  AND  DETRENE  (IE  NECESSARY)  THE  DATA 
CALL  ZEROISE  AN  (CWOPK.NCFM^ 

'*  IN  LOW  THE  DATA 
CALL  FEOLATEICWORK. NOFNT. P> 

CONFUTE  THE  FCURIEP  TRANSFORM 

mofni=ibinafyTncfni  ) 

m  i  i iffispsiiw1' 

OUTPUT  RESULTS 
rp=(rcFNT/2T 

CALL, REALP ART (EAT A .CWCRK .MCPNT 

i?=!5 

FREC=X 1/MOPNT  v 

IF  CIATAII+l'. GT.-40. )  CC  TC  45 

lataTi+i)=-4C. 

*3111(16.46)  FHEC.IATA  HI) 
IOBKATtlX,G15.7,10X,G15.7* 

‘ “INLE 


45 
n 

END 

SUFROUT IN  E  RANGE  (X  ,  NS  ,  XN  AX  ,XMK  .IX) 

PROGRAM  DETERMINES  MAX, MIN  VALUES  AND  ANPI.ITUEF  Cl  A 

wmm 

XNAX=XMIN 


finjiMmi)  Emu] 

50  CONTINUE 
DX=XMAX-XMIN 

IHM* 

SUBROUTINE  MEAN  ( X  .NS  .AVG  'I 
SUBROUTINE  DETERMINES  MEAN 
IIMENSION  X (500; 

SUM=0 . 0 


VALUE  CE  A  VECTOR 


LEPTH  SCAIE 


,NF,NL) 


VECTOR. 


1182 


Ztii 

tee 

600 

821 

set 

see 

850 

£00 

*01 

1000 

im 

1200 

1210 

1220 

1400 

Ic00 

l£00 

1610 

1  7§^ 

1710 

1766 
1760 
176c 
1756 

1767 

m=. 

1  C  Ct 

1760 

ills 


h  r. 

UNCTIC 


EIPENS  ION  AA(54) 

FABArETEB  SETUP:  V=SEr  BAIE.  XI0=TBACFF  IMPACT  DEPTF 
X1C=40 .0 


WEI?E?6^100 )  D^P.EIP 
ICEt'kU  IEP=  .F10.5.5X,  I 
PI=3. 14159 

¥Se1=i 1=1,64 

il=iKiSianif5»Hr!F"xlc> 

C1=-((XIC+LEP-XI;/V)*(PI*?I* 
D1=“((XIG+DEF“XI  /V  *  4*Fa*P 
E1=-((XI(>IEF-XI)/V  W3*FI*F 
A=EXP(A1J  4 
P  =  Tl  ./LEI )*‘EXF (II J 
C  =  C2 ./DEI JVEXF (Cl  1 
E=  U./CEP  vjx?  El 


CO  TO  20 


r IF-  \F10.5) 


o**  2 ) 

LEF**2) 

DEF**2) 


20  90  21  J=1 EEX , 64 

21  AA ( I;=0 .2 

90  25  1=1,64 
25  XTCT  =XTCI  +  AA ( 1 ) 

1»B  IT  E  (6  >22  )  XTO? 

22  FCEJITI  XTOT-  .Fl?, 
I  c  3  01  1  =  1.64 

2e  AA(l]=AA(l )/XT0T 
PiTUEN 
E.\E 


I 

*22 

eee 

?ee 

eee 

see 

III! 

1322 

1422 

i!ce 

ieee 

i?ee 

leee 

1S02 

22ee 

2122 

22:22 

22ee 

2422 

252e 

2?22 

272k; 

2622 

2Se2 

3222 
3122 
3222 
332e 
34  2  2 


*****#«**************»***«'#*«*«*»?«>?>?« 

iHnbim^PiaittsTtHiittyiunH S!CNAL 


■  a  t  i.  ai  ia  ii  a  i  v  i/s.w  nu«  wwa  A'liiM  •  mu  ■  i  n 

WAVE  IRIQUENCIES  ARE  Bill  AS  (2*PI 

DIMENSION  A (500) ,B(500  )  ,C  (500 \t  (  5 

PI=3. 14159 

AMP=1 . 0 

BMP-1.0 

CMP=1.0 

DMP=0,0 

AE=(2*FI)/21.0  % 

BF=  2*PI  /(8.22 
CI=(?*PI  /(4,62) 


I*  §  4/1  )  A 

UlHil 


GTB):  AF,BE,CF,DF 


1020 


1021 


£f:  lEH  Mis!!’ 

WHITE(6  1000; 

FORMAT i  COMPONENTS  01  THE  INPUT  TEST  SIGNAL  ' ) 

DO  10  1=1 ,NS 
A(I]=AMP*[SIN  AE*I )) 

WRITE (6,1001 )  AMF.AF 

JgRMATj  WAVE  AMPLITUDE  =  \F6.2,  J  2*PI/WVLGTH  =  \F6.2 


32  Cjl 
WRI 


SIN  BE*I  )) 
WRITE(6,1021)  BNF.BF 
EG  30  1=1, NS  , 
ClIJ=CNF*tsiNfCI*I} ) 
Vn ITE (6 ,1001)  CMF.fcF 


TO  40  1=1, NS  , 

42  D  ( I  =DMP*(SIN (DE*I)) 
WRITeT6,1001)  DMP.DF 
ro  100,1=1, NSX  ,  %  ,  % 

100  l(I}»A(I)+fi(I)+C(I)+D(I) 
RETURN 
LNI 


!ra45Kllp!i?S%.*«*N) 

LINENS  ION  A ( 1 ) 

AMEAN=0.0 


EC, 10  1=1, N 
10  jMEM=AMEAj;|+A(I  , 

20  MI j=A ( I  ;-AMEA.N 
RETURN 
END 


600 

900 

iee0 
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